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Dynamical tunneling between symmetry related invariant tori is studied in the near-integrable 
regime. Using the kicked Harper model as an illustration, we show that the exponential decay 
of the wave functions in the classically forbidden region is modified due to coupling processes that 
are mediated by classical resonances. This mechanism leads to a substantial deviation of the split- 
ting between quasi-degenerate eigenvalues from the purely exponential decrease with 1/h obtained 
for the integrable system. A simple semiclassical framework, which takes into account the effect of 
the resonance substructure on the KAM tori, allows to quantitatively reproduce the behavior of the 
eigenvalue splittings. 



Despite its genuinely quantal character, tunneling is 
strongly influenced, if not entirely governed, by the struc- 
ture of the underlying classical phase space 0. Chang- 
ing the dynamics from integrable to chaotic substantially 
modifies the tunnel coupling between two wells, not only 
if the chaos affects the classical motion within the wells 
(as considered in 0), but also if the destruction of in- 
variant tori is entirely restricted to phase space domains 
far inside the classically forbidden region. 

Consider, for instance, a classical system exhibiting, 
due to a discrete symmetry, two congruent but separate 
regular regions. Semiclassical EBK quantization enables 
one to construct "quasi-modes" within each of these re- 
gions (i.e., wave functions fulfilling the Schrodinger equa- 
tion to all powers of %), with exactly the same energy for 
any pair of symmetry related quantized tori. This de- 
generacy becomes eventually lifted by the non-classical 
coupling between the quasi-modes, which selects the sym- 
metric and the antisymmetric linear combination of these 
states as "true" eigenstates of the quantum system. 

This "dynamical tunneling" process Q arises in both 
integrable and non- integrable systems, but with substan- 
tially different effectiveness. Introducing an appreciable 
chaotic layer in between the two regular regions signifi- 
cantly enhances the tunnel coupling as compared to the 
integrable case and induces a huge sensitivity of the cou- 
pling with respect to variations of external parameters 
|^-^|. This phenomenon was successfully explained by 
the interaction of the regular quasi-modes with quantum 
states living semiclassically within the chaotic domain, 
which, due to their delocalized nature, assist at the con- 
nection between these quasi-modes |^~|^]. Based on this 
picture, random matrix descriptions of the chaotic part 
of phase space were shown to reproduce the statistical 
properties of the tunneling rates [||-^) . Recent wave chaos 
experiments on optical and microwave cavities ]Tl| ] 
confirm the relevance of chaos for the tunnel coupling 
in nonintegrable systems. For the classically forbidden 
component of the underlying coupling process, however, 
namely the continuation of the wave function from the 



quantized torus to the chaos border, a general semiclas- 
sical theory is still lacking. 

In the present contribution, we shall not tackle this 
problem directly for strongly mixed dynamics, but con- 
sider a simpler situation, namely the nearly integrable 
regime for which chaos is not appreciably developed. Yet, 
we shall see that, despite a seemingly "regular" phase 
space, the tunnel decay of the wave function and the 
associated coupling rates are non-trivial and cannot be 
reproduced by an integrable approximation of the dy- 
namics. To explain this, we shall emphasize the role of 
classical resonances which come into play through two 
important, but distinct aspects. The first one, the im- 
portance of which in this context has already been recog- 
nized by Bonci and coworkers [Q , is the existence of level 
crossings induced by a resonant torus in the semiclassical 
spectrum. Alone, however, they do not induce any mod- 
ifications, as they arise also in integrable systems. The 
second aspect concerns the coupling of near-degenerate 
levels, which is related to the influence of the destruc- 
tion of resonant tori on the nearby phase space p3| . [l4| . 
The main purpose of this letter is to demonstrate that 
the combination of these two aspects form the basis of a 
mechanism which, via the exponentially decaying tail of 
the eigenstates, controls the tunneling between symmet- 
ric islands. At rather large values of h, this process may 
imply only one resonance. Deeper in the semiclassical 
regime, however, several of them are generally involved. 

To make things more concrete, we shall illustrate our 
discussion on a specific, one-dimensional time periodic, 
example, namely the kicked Harper [15) . Classically, this 
system is described by the Hamiltonian 
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r 5(t — nr) cos q 
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where r represents the kick period as well as the per- 
turbation strength. Fig. |j](al) shows the near-integrable 
classical phase space of the corresponding map for t = 1 
within the fundamental domain — n < q,p < n. 

The quantization of the Kicked Harper map is conve- 
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FIG. 1. (al) Classical phase space of the Kicked Harper at 
t = 1 . Boldly marked are the quantized torus with the action 
Iq = w/6 as well as the island chain corresponding to the 
10:1 resonance, (bl) Phase space generated by the integrable 
approximation H (up to C(r 6 )) of the Kicked Harper map. 
(a2, a3) Probability densities of Kicked Harper eigenfunctions 
at t = 1: (a2) 5th excited state for TV = 2-rr/h = 66; (a3) 10th 
excited state for TV = 126. (b2, b3) Eigenfunctions of H 
at r = 1, weighted by their overlap with the Kicked Harper 
eigenstate in (a2) and (a3), respectively: (b2) 5th (solid line) 
and 15th excited state (dotted line) for TV = 66; (b3) 10th 
(solid line), 26th (upper dotted line), 36th (middle dotted 
line), and 50th excited state (lower dotted line) for TV = 126. 



niently described by the propagator U over one kick pe- 
riod, which is given by the product of the two unitary 
operators, U = exp (—irft' 1 cos pj exp (—irfi~ cos qj , 
representing the kick and the free propagation between 
the kicks, respectively. For ft — 2n/N with integer TV, 
the 2tt periodicity in q and p permits writing the eigen- 
functions ip as Bloch-like functions in both position and 
momentum - e.g., by imposing ip(q + 2ir) = ±ip(q) and 
tp(j)+2ir) = tp(p) with ip the Fourier transform of tp. This 
effectively reduces the eigenvalue problem Uip = e i(p ip 
to the fundamental domain, yielding, for each particu- 
lar choice of the periodicity conditions, a finite spectrum 
with TV discrete eigenphases tp n . 

Fig. ||(a2) shows the probability density of a typical 
eigenfunction in the near-integrable regime at r = 1, 
for ft — 2ir/N with TV = 66. The wave function cor- 
responds to the 5th excited state within the central reg- 
ular region around (q,p) — (0,0) and is localized on the 
classical Kolmogorov-Arnold-Moser (KAM) torus with 
the action Iq = ir/6 (boldly marked in Fig. [l]). We 
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FIG. 2. As a function of TV = 2n/h, we plot the difference 
between the eigenphases of the periodic and the anti-periodic 
wave function [defined by t[>(q + 2n) — ±tfj(q)] for the nth ex- 
cited state at TV = 6(2n + l). Solid line: exact quantum result 
for the Kicked Harper. Dashed line: eigenphase splittings of 
the corresponding states of the integrable approximation H . 
Dotted line: semiclassical result obtained by taking into ac- 
count the effect of the 8:1, the 10:1, and the 14:1 resonance. 



see that the eigenfunction exhibits the usual oscillatory 
structure within the torus and decreases exponentially 
beyond the caustics. However, this tunneling decay is 
not monotonous, but interrupted by a local "shoulder" 
at q ~ ±0.77r. For smaller ft the same description applies 
except that several shoulders usually appear. See, e.g., 
the 10th excited state for TV = 126 in Fig. 0(a3). 

The modification in the tunneling tail of the eigenfunc- 
tions has appreciable consequences for the eigenphase 
splitting Aip n = | tpn — (fin | between the "periodic" and 
the "anti-periodic" state [defined by -0(g + 27r) = ±ip(q)], 
which is, in analogy to the energy difference between 
symmetric and antisymmetric states in double well prob- 
lems, entirely dominated by the tunneling tails of the 
wave function. Fig. || shows, as a function of TV = 2ir/ft, 
the splittings for the nth excited states at TV = 6(2n + 1), 
which are all localized on the torus with action Iq . A clear 
departure from a purely exponential decay with 1/ft, to 
be expected for integrable dynamics, is observed. 

In order to reveal the cause of this perturbation, we 
expand our wave function in terms of eigenstates of the 
time-independent Hamiltonian H that is best adapted to 
the KAM structure of the regular region. This effective 
Hamiltonian can generally be expressed as an asymptotic 
series in the perturbation strength r. We use, in our case, 
a method introduced by Sokolov |L6| and obtain 

H = cosp + cos q — — sinp sing + 0(t 2 ) (2) 

in first order in r. Fig. [l](bl) shows the phase space 
generated by the integrable approximation H up to the 
6th order in r. We see that the KAM tori of the Kicked 
Harper are well reproduced within the regular region. 
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Calculating the overlap \(4>m\4'n} \ between the Kicked 
Harper eigenstate ip n and the mth excited eigenstate of 
H identifies the structure in the tunneling tail in terms 
of contributions from the eigenstates of the integrable 
approximation. In fact, this overlap is, as expected, al- 
most unity at to = n, decreases rapidly with \m — n\, and 
exhibits relatively pronounced local maxima for particu- 
lar highly excited states of H . Plotting their probability 
densities [dotted lines in Fig. |T](b2,b3)] clearly confirms 
that it is the admixture of these particular components 
which modifies the tunneling decay of ip n . 

This "selection rule" suggests a resonance phenomenon 
as underlying mechanism, principal aspects of which have 
already been discussed in the literature (l^Jl^l . The clas- 
sical phase space of the Kicked Harper map exhibits is- 
land chains corresponding to nonlinear resonances be- 
tween the kick and the unperturbed dynamics, two of 
which, the 10:1 and the 14:1 resonance, are visible in 
Fig. 0(al). A r:s resonance, with r islands, generates in 
its neighborhood a periodic perturbation of H and thereby 
couples ip n to ipn+rk with integer k, and it is precisely 
the combination of these couplings that gives rise to the 
shoulder structures in the tunneling tail. 

To demonstrate that the above interpretation is basi- 
cally correct, we use a semiclassical description of these 
couplings between quasi-states to reconstruct the wave 
function in the tunneling regime and to reproduce the 
behavior of the fluctuations in the eigenphase splittings. 
Following the lines of the derivation undertaken by Ozo- 
rio de Almeida |Q, we perform our semiclassical study 
in the action-angle variable representation (/, 9) associ- 
ated with the integrable dynamics of H. In this repre- 
sentation, the Kicked Harper Hamiltonian is written as 
H = H (I) + V(I, 9, t) with H Q {I) = H(p, q) the unper- 
turbed integrable Hamiltonian. 

Let us, for the moment, assume that the phase space 
of the integrable system is perturbed by only one single 
island chain corresponding to a r:s resonance. Follow- 
ing standard secular perturbation theory |]l7| , we place 
ourselves in the co-rotating frame by the transformation 
9^9 = 9 — lnstjrT to the slowly varying angle 9. It 
is then possible to perform a canonical transformation 
towards new action- angle variables (/*, 9*) for which the 
perturbation V is effectively replaced by its time-average 
% s = (1/rr) J Q rT dtV(I*,6* + 2irst/(rT),t) over r kick 
periods. This yields a time-independent Hamiltonian 



H (I*) 



2ns 



I* + V r:s (I*,6*), 



(3) 



for the effective dynamics in the vicinity of the r:s res- 
onance. The potential V r:s is given by a 2ir/r periodic 
function of 9 and can therefore be written as 



V r:s = 



OO 

E 

fe=i 



V r k :s cos(Ax0* + 



(4) 



In practice, we directly extract the Fourier coefficients 
Vr- S , neglecting their action dependence, from the struc- 
ture of the separatrix associated with the r:s resonance. 
In the action-angle variable representation, the functional 
form litj(9), itp ?(0) of the "upper" and the "lower" sep- 
aratrix between the resonance substructure and the KAM 
domain can be numerically constructed in a straightfor- 
ward way. From the knowledge of Ho(I) and from the 
condition that H r . s (litJ (9), 9) be independent of 9, the 
coefficients Vf s are calculated by the Fourier analysis of 



V r -. s (9) 



H (I) - 



2ttsI 



— const (5) 



J=/ r:3 +A/ s<!p (0) 



with I r - S the action at the r:s resonance and AI sep (9) = 
(I$(0) ~ Isc P \9))/2 [using I r:s + AI scp instead of litj , 
we compensate for the error that is introduced by the 
discrepancy between the (/, 9) and the (/*, 9*) represen- 
tation] . 

The Kicked Harper eigenstates ip n can now be explic- 
itly constructed as linear combinations of the eigenstates 
tp n of the integrable Hamiltonian within the framework 
of the first order perturbation theory, which involves the 
unperturbed eigenenergies E n of ip n and the matrix ele- 
ments {tpn'Y^r-.s^n) of the perturbation V r[s . These ma- 
trix elements are directly calculated in the action-angle 
variable representation (1,9); using Eq. (||) and repre- 
senting the unperturbed eigenfunctions semiclassically as 
(9\tjj n ) — exp(m0)/v27Tj we obtain [a = sgn(n' — n)] 



(lpn'\Vr: S \lpn) = -T^ S \n'~n\ ,kr exp(i (J £fe) . 



(6) 



fe=i 



We notice that nonzero couplings arise only between 
states the quantum numbers of which differ by inte- 
ger multiples of r. As a consequence of this selection 
rule, the perturbative expression for the nth Kicked 
Harper eigenstate involves only contributions ?/v with 
n' — n = kr for integer k, and the eigenphase splittings 
Aip n — \tp n + ^ — fn \ between periodic and anti-periodic 
states are to a very good approximation given by 



Aip n = r™: s fe A£„ +fer with 



V r k s /2 



E n - En+kr + 2nTlsk / T 



(7) 



(and r«f=° = 1). A(p n = \<f { n } - (p n ) \ denotes the 
nth eigenphase splitting for the integrable Hamiltonian 
H and can be semiclassically computed using standard 
WKB techniques for integrable dynamics. In principle, 
the sum on the r.h.s. of (^) involves all couplings to the 
states the quantum numbers of which satisfy the selection 
rule n' — n = kr. In practice, it is typically dominated by 
the contributions with quasi-energies closest to degener- 
acy in the Floquet spectrum of the kicked system. Such 
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a near-degeneracy occurs if the action of the r:s resonant 
torus lies close to the arithmetic mean of the classical 
actions associated with ip n and i/Vi+rfc , in which case the 
energy difference E n — E n+ km is determined by the lin- 
ear approximation of H near the resonant torus, which 
exhibits constant eigenphase spacing 2-ns/r. 

The dotted line in Fig. ^ shows the semiclassically cal- 
culated eigenphase splittings of the states ip n that are 
localized on the torus with action Iq. Below N = 100, 
only the 10:1 resonance is taken into account. Above 
N = 100, multiple coupling schemes involving also the 
8:1 and 14:1 resonances become important. This, how- 
ever, does not pose any essential difficulty. Taking into 
account the effect of each r:s resonance separately, we 
can recursively apply the mechanism described above for 
a single resonance, which amounts to generalize Eq. (Q) 
to a multiple sum over products of I 1 "'*. This sum is 
again dominated by near-degenerate contributions; for 
N = 126 for instance [Fig. ry(a3,b3)], it is given by 

A ^10,2^16,1^36,1 a ~ 

We see in Fig. || that the agreement with the exact 
quantum result is fairly good, though not perfect. Devi- 
ations from the quantum splittings occur due to neglect- 
ing the action dependence of the Fourier coefficients V£ a 
of the potential (Q) and sometimes also due to acciden- 
tal near-degeneracies of ip n with states that are associ- 
ated with the "opposite" regular region centered around 
(q,p) — (tt, it) beyond the separatrix. Nevertheless, the 
main features in the fluctuations of the level splittings 
are well reproduced by our semiclassical description. 

In conclusion, we have given a simple prescription how 
to reproduce the tunneling rate fluctuations in the near- 
integrable regime by the classical structure of a few im- 
portant resonances. The effect of a r:s resonance on the 
tunneling process of a regular eigenstate is twofold. On 
the one hand, the periodic modulation of the KAM torus 
structure due to the resonance induces appreciable cou- 
plings to higher states that are selected by the periodic- 
ity r of the island chain. These couplings become, on the 
other hand, substantially enhanced at near-degeneracies 
in the quasi-energy spectrum. We have explicitly incor- 
porated this mechanism in a perturbative semiclassical 
framework which, through simple iteration, is general- 
ized to multiple coupling schemes involving several reso- 
nances. As underlying picture, we obtain that each r:s 
resonance allows the wave function to "hop" from one 
side of the island chain to the other, with the quantum 
number changing by integer multiples of r. For large fi 
only one or few resonances with relatively small r come 
into play, whereas for small h the tunnel coupling pro- 
ceeds via a succession of hops over several adjacent island 
chains, rather than by one single, but longer, classically 
forbidden event. Comparison with the quantum tunnel- 
ing rates shows good agreement, and confirms in particu- 
lar also that contributions from the chaotic part of phase 
space do not play a role in the near-integrable regime. 



Though exemplified only within the Kicked Harper 
model, we expect that the resonance-assisted coupling 
phenomenon discussed here represents the fundamen- 
tal tunneling mechanism for general near-integrable sys- 
tems. An open and interesting question is to which ex- 
tent the role of the resonances retains its importance 
also in strongly mixed regular-chaotic systems where a 
macroscopic fraction of quantum states is associated with 
chaotic domains. We believe that at least the "regular" 
contribution to chaos-assisted tunneling, i.e., the cou- 
pling from the torus to the chaos border, is now amenable 
to straightforward reproduction via the identification of 
the major resonances inside the regular domain. 
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